Extended hydrodynamics from Enskog's equation for a two-dimensional system 

general formalism 



O 

o 

(N 
> 

o 



B 



I 

O 

o 



> 
oo 
cn 

o 

O 

o 



I 

o 
o 



X 



Hideaki Ugawa and Patricio Cordero 

Departamento de Fi'sica, FCFM, 
Universidad de Chile, Santiago, Chile 
(Dated: February 6, 2008) 

Balance equations are derived from Enskog's kinetic equation for a two-dimensional system of hard disks 
using Grad's moment expansion method. This set of equations constitute an extended hydrodynamics for mod- 
erately dense bi-dimensional fluids. The set of independent hydrodynamic fields in the present formulations 
are: density, velocity, temperature and also — following Grad's original idea — the symmetric and traceless pres- 
sure tensor pij and the heat flux vector q*. An approximation scheme similar in spirit to one made by Grad 
in his original work is made. Once the hydrodynamics is derived it is used to discuss the nature of a simple 
one-dimensional heat conduction problem. It is shown that, not too far from equilibrium, the nonequilibrium 
pressure in this case only depends on the density, temperature and heat flux vector 

PACS numbers: 51.10.+y, 05.20.Jj, 44.10.+i, ()5.70.Ln 



I. INTRODUCTION 

In 1921 Enskog introduced a kinetic theory for dense 
gases lit] which is known to yield a very good approximate 
description of the behavior of gases, particularly transport 
phenomena, as in |01 and references therein. Today En- 
skog's original theory is known as the standard Enskog theory 
(SET) illl 1 because after the pioneer work of van Bei- 
jeren and Ernst Ist] there are several new versions of Enskog's 
theory collectively called revised Enskog's theory (RET) 16]. 
Among the latter there are versions that have been extended 
to describe condensed matter To Navier-Stokes level both 
SET and RET lead to the same results jsl Ht], whether or not 
an external force is present. For SET and RET, using approx- 
imations in section |Vl the same hydrodynamic equations are 
obtained. 



sure. 



Once a kinetic equation is given it is possible to de- 
rive hydrodynamic like equations, thus reducing the degrees 
of freedom from those of the velocity distribution function 
/(r,c,f) to the degrees of freedom of a finite set of hydrody- 
namic fields. Of these methods we mention the Chapman- 
Enskog method lH, |^ IToll and Grad's moment expansion 
method lfTll[l2ll . The first is widely used and illustrated even 
in textbooks. Much less attention has been paid to Grad's 
method which is our central tool in this paper. Both Chapman- 
Enskog and Grad's methods are truncated expansions. 

Perhaps the first effort to obtaining hydrodynamic equa- 
tions (balance equations) from Enskog's theory using Grad's 
moment expansion methods was by Schmidt, et al [13] in 
1981. They did not start from the local distribution function 
but from a Maxwell distribution function as kernel, totally 
linearized their hydrodynamic equations, that is, completely 
eliminated the products of differences between the hydrody- 
namic variables associated to the equilibrium state and the in- 
homogeneous nonequilibrium state. From this they obtained 
approximate hydrodynamic equations such that the conserva- 
tion laws for mass, moment and energy are not exactly satis- 
fied. They also obtained expressions for the transport coeffi- 
cients (viscosity and heat conductivity) and the thermal pres- 



In 1988, Kremer and Rosa fl^l obtained hydrodynamic 
equations from the local equilibrium distribution function as 
kernel linearizing the collision integral in Enskog's equation. 
In particular they introduced an approximation which com- 
pletely eliminates the second order terms of the collision inte- 
gral. They derived sound dispersion relations for monoatomic 
gases using normal mode analysis. 

In 1991, Marques and Kremer ifTsIl developed Kremer and 
Rosa' work lfl4ll and obtained linearized hydrodynamic equa- 
tions involving the second order terms of the collision inte- 
gral. Consequently their equations are less approximate than 
in Ref. 1 1 3il although in the bi-dimensional case, eliminating 
the products between the hydrodynamic variables which are 
deviations from the value at equilibrium, their equations co- 
incide with our results. Furthermore they obtained linearized 
Burnett equations for monoatomic gases where they eliminate 
third order terms from the collision inte gral . 

In 1996, Rangel-Huerta and Velasco lUol obtained the hy- 
drodynamic equations from the local equilibrium distribution 
function as kernel with partial linearization of the collision in- 
tegral in Enskog's^quation. They also eliminated collision 
terms as in Ref. |14f| and obtained the extended equations as 
in Ref. liTsll : in the second cited paper of Ref. 1 1 6il Grad's dis- 
tribution function without approximation is used, whereas in 
Ref. ifTsIl the authors use a Maxwellian distribution function 
to get the second order collision term. Therefore, the terms 
of the second order spatial derivative of the pressure tensor 
and heat flow vector are added in the balance equations for 
the pressure tensor and heat flow vector. They derived gen- 
eralized transport coefficients wj^th the methods used widely 
in generalized hydrodynamics I . In 1997, Ref. liTsIl devel- 
ops the extension of Grad's moment method for the RET, and 
in 1998, obtained hydrodynamic equations which are applica- 
ble only to a system undergoing uniform shear flow lfl9ll . In 
2001, Rangel-Huerta and Velasco ll20ll eliminated the second 
order terms of the spatial derivative of the pressure tensor and 
heat flow vector from the balance equations for the pressure 
tensor and heat flow vector of the former results lfl6ll and ob- 
tained a similar but extended hydrodynamic equations similar 
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to Ref. ifTsIl nevertheless, get better agreement with molecular 
dynamic simulations ll2lll . 

The present article provides extended hydrodynamic equa- 
tions derived from Enskog's equation using Grad's moment 
expansion method in the bi-dimensional case. They are more 
complete than a linear approximation but still they are the re- 
sult of an approximation scheme that we later explain. As far 
as we know this is the first time that Grad's method used to ob- 
tain extended hydrodynmaic equations from Enskog's equa- 
tion has been published beyond the a linear approximation in 
two or three dimensions. 

In addition, we will apply this hydrodynamic equations to 
a one dimensional steady state heat conduction case. There 
are many studies on this theme: for example, the experi- 
mental investigation in Ref. ll22ll and the theoretical one in 
Refs. 1231 124I1 . Gross and Ziering ll25ll investigated it start- 
ing from Boltzmann's equation with Grad's moment method. 
To our knowledge, no one has investigated it starting from 
Enskog's equation using Grad's moment method. Recently 
Kim and Hayakawa 12611 investigated it starting from Boltz- 
mann's equation using the Chapman-Enskog method for hard 
core particles, Maxwellian ones and BGK's ones ll27ll and 
Hayakawa et al ll28ll investigated it for hard disks using 
molecular-dynamic simulations. Furthermore, they tried a 
test of the nonequilibrium steady state thermodynamics (SST) 
proposed by Sasa and Tasaki 12911 and criticized it. The state 
defined in Ref.|f25] for a gas in a nonequilibrium steady state 
is a one-dimensional heat conductive configuration in contact 
with an equilibrium state through a special porous wall (called 
the perfect jx wall ll29ll ). There is a nontrivial pressure differ- 
ence between the equilibrium and nonequilibrium state paral- 
lel to the direction of the heat flow. We have considered this 
phenomenon with our hydrodynamic equations in Ref. Isoll . 

Parenthetically we mention that it is meaningful to check, 
using kinetic theory, whether one of the main assumptions 
of SST ll29tl — the nonequilibrium steady state pressure — is a 
functional of the density, temperature and heat flow. We find 
that, in fact, it is correct provided the system is not too far 
from equilibrium. This is consistent with the range predicted 
in Ref. il29ll for this nontrivial pressure difference to appear 
Note that one of the useful merits of Grad's moment method 
is that one can easily find relations between nonequilibrium 
thermodynamic variables. 

The organization of this paper is as follows. In Sec. HH 
Grad's moment expansion method is introduced. Sectionlllllis 
about Enskog's equation. In Sec.|IV]we formulate the macro- 
scopic balance equations. In Sec. |V] the simplified collision 
contribution used in this paper is constructed showing that, 
for the approximations introduced here, there is no difference 
between the macroscopic balance equations obtained by RET 
and SET. In SeclVIIthe hydrodynamic equations that our pro- 
cedure yields are presented and the collision frequency coeffi- 
cient X used in this paper are introduced. In Sec. lVIII we apply 
our hydrodynamic equations to a simple heat conduction case 
in which there is a steady heat current between two parallel 
plates at slightly different temperatures. The closed expres- 
sion of the pressure only depends on the hydrodynamic fields: 
density, temperature and heat vector. A final remarks on this 



example are given. In Sec. lVIIII we summarize our results. 



II. THE VELOCITY DISTRIBUTION FUNCTION AND ITS 
EXPANSION 

To define dimensionless variables A from the physical vari- 
ables A we choose the mass of the particles to be the unit of 
mass, and the temperature is measured in energy units so that 
Boltzmann's constant is 1 . Further, we take a reference dis- 
tance L and a reference temperature Tq so that the dimension- 
less coordinates x, velocities c and time are such that 



t = 



'-x/L, c = c/v7b 



To, 



L 



Assuming that the system is in a rectangular box of size 
X Ly the hydrodynamic fields can be rescaled defining di- 
mensionless fields as follows. 



V = 



q 

F 



-n 

— T 

To 

LyLv — 

- p 



NTr 



To 



3/2 



number density, 
hydrodynamic velocity, 
temperature, 
pressure tensor, 
energy flux, 

external force, 



(1) 



The distribution function / is replaced by a dimensionless dis- 
tribution function /, 

■' ~ ■'' 

In the following sections some auxiliary dimensionless 
quantities will be used 



a,. 
Po 

Kn 

5 



L,Ly 4 

8 V2 ^ 
- =Knpo 



(2) 



Knudsen number, 



where a is the disk's diameter, is the number of disks and 
£ = is the mean free path at equilibrium. In these units, 

for example, the free flight time for Boltzmann's gases at equi- 



librium temperature Tq is 



/ JiKn 



L is any macroscopic distance 



that we take as relevant depending on the particular problem 
to be analyzed. For the purposes of the present paper we will 
typically assume that L is of the order of Ly. 
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Following Grad, we make a formal eight moments expan- 
sion of the distribution function ifTlll 



fir,c,t) = /m^g, 



where 



/m = n(r,f) 



27tT{r,t) 



exp 



2T{v,t) 



(3) 



(4) 



In Grad's method the factor is written as an expansion in 
Hermite polynomials of the peculiar velocity C = c — v(r, t). 



4>a=E«r(r,0<'(C), 

n=0 



(5) 



where c is the molecular velocity and where H^"^ is a tensor 
with n subscripts, / = {i\ ■ ■ •/„) as well a polynomial of n-th 

in) 

degree in the components of C. Also the coefficients a] are 
tensors of order n. The eight moments expansion corresponds 
to truncating the expansion in Eq.(|5]) up to n ~ 3. In dimen- 
sion 3 this leads be the well known 13 moments methods used 
by Grad himself working upong Boltzmann's equation. It is 
widely used. 

Grad's method consists of first replacing / given by Eq.(l3]i 
but with the truncated expresion for Og, in the kinetic equa- 
tion (Boltzmann's or Enskog's), and then projecting the ki- 
netic equation to the different Hermite polynomials used in 
the tructated expansion. This procedure — deviced by Grad — 
leads to obtaining as many balance equations as polynomials 
were included in the truncated expansion for Og. 

The Hermite polynomials satisfy the orthogonality relations 



where i,] =1,2. Instead of using the full set of third order 
Hermite polynomial Hf-j^ we use — as Grad did — the corre- 
sponding contracted Hermite polynomials Hf defined by 

Ci (c'- 



W 



Vf\T 



(11) 



The coefficients a" in Eq.© can be expressed in terms of 
the hydrodynamic fields n(r,f), T{r,t), and the flux vectors 
associated with the velocity and the energy, called the kinetic 
part of pressure tensor and the energy flux: i^^(r,f), q^{r,t) 
according to the following sum rules, 



fdc = n(r,f), 
C(r,t)/dc = 0, 



ic(r,t)V«'c = «(r,t)r(r,f), 



J Ci{r,t)Cjir,t)fdc = 



/ 



C{r,t)C^{r,t)fdc = qHr,f), 



(12) 
(13) 
(14) 
(15) 
(16) 



where pjj is the traceless part of Pj'j, 

Ff.^n{r,t)T{v,t)5,j+pij{v,t). 



The factor <I>g for the eight moments method turns out to 



be 



1 



2kT 



2T 



dC = n\5„,n8. 



and they are, for example, 

i/« = 1, 
1 _ Ci 



H. 



CiCj 



QC jCk {CiSjk + Cj5ki + Ckdij) 



Ti 



(6) 

(7) 
(8) 

(9) 

(10) 



*r7 = l 



2«r2 



1 

PijCiCj + -^^ 



4T 



C q' (17) 



We refer to the "kinetic part" of the fluxes to stress that 
kinetic theory, in principle, includes contributions to the fluxes 
associated to the interaction between the particles. 



m. ENSKOG'S EQUATION 



Enskog's equation is 



+ci • Vi +F • VeJ/(ri,ci,r) - ^ / [Z(ri,ri + 5k |«)/(ri,ci,f)/(ri + 5k,4,f) 

-X (ri, ri - 5k I «)/(ri , ci , f )/(ri - 5k, C2 , f )] (g ' k) 0k dkdcj (18) 
I 



where k is the unit vector from the disk centered at ri + 5k 
to the disk with center at ri upon collision and is is inte- 
grated over all the unit vectors while g = C2 — Ci and 0^ is 



the Heaviside function 9^ = 9{g-k) and where the two ve- 
locities of the disks after collision are c\ = Ci + (g • k)k and 
c'^ = C2 — (g • k)k, respectively. 
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The collision frequency x(ri,ri ± 5k |«) is the pair distri- 
bution function of two hard disks at contact and n is the num- 
ber density defined by Eq. (fT2] i . This x is given by the proce- 
dure of statistical mechanics, that is, by the virial expansions 



of the pair distribution function and of course its expression 
depends on whether one uses SET or RET. 

Concretely, with r' = r + 5k, 



J 



r^^(r,r'|«) 



l+«(i(r + r')) J V(r|r3)t/r3 + ^«2(l(r + r')) J y(r,r'|r3r4Mr3t/r4 - 

3 3^ 
= Xc + 2 k • Vz, + y kk : VVXc + ■■■ 

;t*^^(r,r'|n) = 1+ y"«(r3)y (r,r'|r3)<ir3 + ^ y"«(r3)«(r4)y (r,r'|r3r4Mr3^/r4 + • • • 

;t, + |^(k.V„) + ^^(k.V„)(k.V«) + ... 
8 on^ 



(19) 



(20) 



2 dn 



where 



Xc=l+n{r) / y(r,r'|r3Vr3- 



2! 
I — 



y(r,r'|r3r4)t/r3c/r4- 



(21) 



and where the usual Hushimi V- function is introduced Isi lsTll . 

Expanding the right hand side inpowers of 5 up to second 
order following the Appendix of |8y. 



rhs 
rhs 



SET 



SET 



■^0 

Jo 



SET ,_gjSET^g2jSET^ 



RET^gjRET^gljRET 



where 



and 
^0 



ikxcSifJ^ - /1/2) (g ■ k)ekrfkrfc2, 



^11 = ikxc!^- (/i'v/^+/iV/2)(g-k)ekdkdc2, 

= ^/(k-V;f,)(/i'/2+/i/2)(8-k)MkdC2, 



-Xc I kk : (/; VV/^ - /i VV/2 ) (g • k) ekrfkdcz , 



J22 = ^/(k-V;t,)k(/;V/^-/iV/2)(g-k)ekdkdc2, 
J23 = 25k/(kk : VV;t,)(/i72 -/i/2)(g-k)MkdC2, 



J24 



^^J{kk: VVn(r, t) ) (/[/^ - .A /2) (g • k) Ot^dkdc^ 

(25) 

where the following abbreviations; fi ~ /(ri,Ci,f), /[ = 
/(ri,ci,f), /2 =/(ri,C2,f) and f!, = /(ri,4,f) are intro- 



(22) 
(23) 



jSET 


jRET 

— Jo 






tSET 


jRET 
— Jl 


= Jn +-/12, 




tSET 

•ll 


= J2I + 


/22 +-/23, 




jRET 
J2 


jSET 
— ■'2 


+ J24, 


(24) 



duced. Here Xc is concretely given by the virial expan- 
sions Jit], but, as seen below, we do not use such an expression. 
IV. THE BALANCE EQUATIONS 

The macroscopic balance equations are derived following 
Grad's prescription, successively projecting the kinetic equa- 
tion to the different Hermite polynomials being used. 

If Enskog's equation is multiplied by a function 0i of the 
peculiar velocity Ci = Ci — v and integrating over Ci yields 



0i^/ic/ci + J (j)iCi-Vfidci + j (piF -V.Jidci 



= 1 



tSET 



tRET 



(for set) 
(for ret) 



where 



/,■ = j (PiJidci (!= 0,1,2). 



(26) 



(27) 



The balance equations of our interest, in particular the ba- 
sic conservation equations (mass, momentum and energy), are 
obtained multiplying Enskog's equation by = l,Ci, jC^. In 
such cases 0i + 02 = ^i' + 02 i^ (prime indicates post col- 
lision velocities) so that 



tSET 



tRET 



0, 



(28) 



and 
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tSET _ jRET 

h - h 



TrKn 



V • j Xci^i - 0i')/i/2(g • k)9kkdkdcidc2, 



TrKn 



k-/i/2Vlog^ 

/2. 



(g • k)0kkc/kc/cit/C2. 



In the case when the external force is uniform the general conservation equation may be written as 

d 



dt 



(/)i/i£/ci + V. / (/)iCi/i£/ci-F- / 0iV,,/ifl!ci = -V.^ 



where /= - V -^. 

In the case when is not associated to a microscopic conservation law the results are 

jSET ^ jRET = i|- I (,/,( - </,i)/l/2(g • k)(kdkdCidC2, 



jSET „ jRET 



tSET 



'1 



TrKn 



'1 



Xck-flVf2 + -ik-V,Xc)flf2 



(k- V;t,)k(/i V/2) +x,kk : /i VV/2 + -(kk : VVx,)/i/2 



(g • k)9kdkdcidc2, 

(gk)0k(ikiici^fC2, 



rRET 



TrKn 



(k • Vx,)k(/i V/2 ) + ;t,kk : /i VV/2 + 



1 dhc 
4 5«2 



(kV«(r,0)(kV«(r,0)/i/2 



r 



(29) 
(30) 

(31) 

(32) 
(33) 
(34) 



(g-k)0k'ikc/cit/C2. (35) 



V. SIMPLIFIED COLLISION CONTRIBUTION 

The collisional term in both Boltzmann's and Enskog's 
equations is quadratic in the distribution function, hence there 
are integrals which involve products H'" which have to be 
projected to different //"s. Even in the case of the 13 moments 
method for Boltzmann's equation with the contraction as in 
Eq. ( fTTI) . Grad proposed to consider only those contributions 
such that ^ + m = « even though there is no clear mathematical 
criterion to make such simplification except for Maxwellian 
molecules because this case the collision term is integrable 
and such quadratic terms do not appear 

Nevertheless, Grad was able to prove that very little error 
is made even in the case of non-Maxwellian molecules by re- 
stricting the contributions to £ + m = n ifTTl \\% . Although we 
do not have such a prove in the case of Ebskog's equation we 
follow the same prescription. 

In our present case, to be able to obtain general hydrody- 
namic equations — following the present procedure — which 
are not overwhelmingly complex it is necessary to make some 
further simplifications. 

Since <I>c is a correcting factor we may consider that 



its contribution can be regarded as 1 plus something small, 
schematically <I>g = 1 + e. Since the colHsion integral is 
quadratic in the distribution function, the integrand has a prod- 
uct <t>(3i<t>G2 = 1 + El + £2 + eie2- Our simplificatory strategy 
will be first that: to evaluate Jq we keep all the contributions; 
to evaluate J\ we keep the "order e" terms, namely, <t>Gi'J'G2 ~ 
<t>Gi +^G2 — 1; to evaluate J2 we take (^g\^g2 ~ 1- 

Therefore, in Eq. ( |33] |, the following relations are required, 

/1/2 =/mi/m2(*G1+*G2-1), 

/l V/2 = /mi (V/m2)(*G1 +*G2 - 1) +/mi/m2V<I>g2 

and in Eqs. (l34b and (1351 ) we use 



/1/2 = /mi/m2- 



(36) 



For hard disks, /mi/m2 = fmfm- Hence, 724 in Eq. ^ van- 
ishes, and from Eq. (|24] |. the SET's collision term coincides 
with the RET's. Additionally, also ^23 in Eq. (l25b vanishes for 
the same reason. Consequently, Eq. ( |35] ) is reduced ( |34| ) and 



jSET ^ jRET = - 01 ) [(k • V;t,)k(/Mi V/m2) +Xckk : /mi VV/M2] (g •k)Mk«'ciflfC2. (37) 



After introducing the approximations mentioned above, we directly calculate, using Maple routines, the collisional term 
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in Enskog's equation for all contributions originated in the 
products H'H"\ 



VI. THE HYDRODYNAMIC EQUATIONS 



and the shear tensor S which is defined by 

'■' dxi dxj dxi 



Sij, 



(46) 



The central result of the present paper are the hydrodynamic 
equations that we have derived from Enskog's equation ac- 
cording to what has been explained in the previous sections. 

To the best of our knowledge, these hydrodynamic equa- 
tions are published for the first time. For the two dimensional 
case we have not encountered any paper at all while for three 
dimensional case the published literature refers to linearized 
versions. For the sake of clarity the balance equations for con- 
served quantities (mass, momentum and energy) are written in 
terms of the hydrodynamic quantities: density, hydrodynamic 
velocity, temperature, pressure tensor and heat flux; whereas 
the balance equations associated to the pressure tensor and the 
heat flux are written in terms of the symmetric pressure tensor 
and heat flux vector. This is what is also done in almost all 
well-known articles for the extremely dilute system. 

We use the following notation as Euler's (material) deriva- 
tive; ^ = ^ + ^ ' Summation is implied with repeated 
tensorial indeces. 

The final extended hydrodynamic equations are the follow- 



ing: 



Dn 

— = -nV-v 
Dt 



(38) 



where Q = j3fyq'' + qcol, 



Bo = 1 + S Ten, 



qcol = -S'-^xy^'^T (47) 



and where the kinetic contribution to the heat flux vector Q is given 
by q*^ following the definition Eq. iT2\ . Additionally, the remain- 
ing terms, which give the collisional contribution to the extended 
hydrodynamic equations, are given as below. Ip is split in the form 
Ip = Ip + Ip+Ip where 



0rKn 



XcVf 



"Pij + j^[2q'lq'j-^UlSij\ 



s_ 

Kn 



u dn u dn u dn 



4 



d 

dxi 



' dxj dxi 
d , _d_ 

OXl 



k^Xc ^ k dXc k^Xc^' 
dxj ' dx; dx\ 



1v, 



-Xc^Pij 

OXl 



;Xcn'TS: 



(48) 



(49) 



D\ 

n — =/iF- V-P 
Dt 



DT 

n = -(V-Q + P: Vv) 

Dt ^ ' 



(39) 
(40) 



D 
Dt 



Pij 



dvi 
h 
d 



dvi 

OXl axi OXl 

1 
2 



dx 
d 



dv,„ 
Plm^ — Oil 
dxi 



dx, 



-nTSi, 



= In (41) 



D 
Dt^ 



■1i 



dx'i 
dx, 



d 

axj 
dT 



1 dvi 
1 d 



-1l 



'2nT 



dxi 



OXl 

d T 

-2^Pli--^PmlPli 

axi n ax„ 

In 



where the pressure tensor P is given as 

P = j3;,p + Peol, 

with 

5 



Pcol 



1 



Kn 



Xcn, 



(43) 



(44) 



Kn 



0rKn 



Xcn^Vf[^ + 2{V-y)\\A5) 



and the 8^ contribution is 



Kn 



2XcnVT 

+ Xcn^-^ 
+ n^s/T 



dn dT dn dT dn dT ^ 
dx, dxj dxj dxi dxj dxj 

'^dT dT dn dT ^ 
dxi dxj dxi dxi ''' 



d^T 



dxidxj dxidxj 



dT dxc ^ dT dxc dT dXcg ' 

dxj dxj dxj dxj dxi dxi '■' 



+ 2Xcn-^^f^S, 



Iq is split Iq = Iq + Iq + Iq and 

I« = -^^[4nrgf 
' %/7rfKn L ' 



(50) 



(51) 



(42) lJ = i|[(l)5 + (2)5 + (3)5] 



(1) 5 

(2) 5 

(3) 5 



d dn 
IXc-nT —pii + ITXc^Pli 
OX) clxi 

, n ^^Xc , , dT 
^■JnT^p,i+4Xc^Pli 

OXl OXl 



^XcnT' 



dxi 



-TXcn'T 



.dT_ 
dxi 



-4n'-T^ 



dXc 
dxi 



:Xc'i 



(52) 
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and where 



.dT 



dxi 



dxi 



dn df^ 
- 36XcTj—Sii + ISXc'iTj-^Vi 
axi axi dxi 



+ ?>6XcnT 



d dvi 
dxj dxi 



(53) 



The above equations have the well known form which comes from 
very general arguments from continuous fluid matter and where es- 
tablished in the XIX century. The detail structure of some of the 
terms comes from kinetic theory. 

Of course it is interesting to compare our results to the ones ob- 
tained by means of the Chapman-Enskog method. For example, Kim 
and Hayakawa (2^ have derived the explicit velocity distribution 
function of Bolzmann's equation using Chapman-Enskog expansions 
to Burnett's level, although restricted to steady states. Also the lin- 
earized Burnett equations in the SET (in the first Enskog approx- 
imation) have been used by Alves and &emer (32ll to study light 
scattering from density fluctuations. However it would be difficult 
to compare our approximations with those from Chapman-Enskog's 
method. They are quite different methods indeed. 

Furthermore, in order to calculate any physical quantity, we first 
determine Xc- The static pressure p is given as p = jTrP. 

The static pressure p is given as p = jTrP. This automatically 
yields that the equation of state is given by Eqs. l |43t and i l44b . Then, 
the relation between Xc and the static pressure p is given by 



Xc 



p/nT-l 
2P() 



(54) 



where the relation; 5 = poKn is used. In general, Xc is related to 
the pair distribution function Q]. In the case without external force, 
the equation of state at equilibrium is very well approximated using 
Henderson's expression 13311 



_P_ 

nT 



P6 



(1-PoP 



(55) 



which is what we use. 

In the case with F ^ 0, SET does not yield the correct single- 
particle equilibrium distribution function, whereas RET do (34ll . 



VII. HEAT CONDUCTION BETWEEN PARALLEL 
PLATES UNDER A STEADY STATE CONDITION 

A. Reduced hydrodynamics 

The above hydrodynamics should in principle allow us to study 
a wide variety of dynamic problems involving dense gases. In this 
section we look at a quite simple situation in this wide context. It 
is the case of a one dimensional heat conduction case between two 
parallel plates. 

The system consists of gas between two infinite parallel plates 1 
and 2 separated by a distance L. The plates have fixed temperatures 
Ti and T2 respectively. A schematic representation in Fig. [T] shows 
the Y axis defined perpendicular to the plates while the X axis is 
placed on plate 1 . 

In general, under circumstances with a gradient of temperature 
perpendicular to the walls, there is a difference between the temper- 
ature of the gas by the plate and the temperature of the plate itself. 





Y 






T2 








Ti 




\- 


X 



FIG. I: Configuration of the simple heat conduction system as de- 
scribed in the text. 



This is a well-known effect called thermal slip. For simplicity's sake 
we neglect such difference. 

For this simple system the extended hydrodynanric equations dras- 
tically simplify. The basic concrete equations solved here are the 
following: 

Pyy{y) = Pyy = COnSt., 2v(.>') = Qy = COUSt., 

Pxy{y)=Pxy{y)=Q, 



1 dq\{y) 

2 dy 



Pyy - 

Qy = 
8 

0fKn 

JL 

i/7rKn 



Pxv (y) + 



,3 ,3 



^.[y)-8' Xcn{yY^T{y) 



dT(y) 



'Ay)pxx{y)- 



./TTKn 

q'yiyf 



128r(y)^ 



8 

+ 4Kn^^ 



dy 



5g (y) 



Xcn{y)^/T^ 



, dn{y)dT(y) 1 ^Cy) f dn{y) 
dy dy 2 T{y) \ dy 



+«(y) 



dy 
d^T(y) 



My)- 



dy 



dy 



(56) 
(57) 



(58) 



-T{y) 



dpxxjy) ^ dp_,,y{y) dT{y) 
dy dy dy 



+ 2n{y)T{y) 



dpxxiy) , T{y) , sdT(y) Pxx{y) dp^xiy) 



+ —r-rPxxiy) 



dY n(y) 



dy i{y) dy 



^J7tT{y)}^ 



Hy)T{y)-pAyH{y) 
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2 dy 



7 , , , , dn(y) , , dT(y) 



dy 



--ln{YYT(y) 



dT{y) 



dy 



(59) 



r 



B. Functional expression of the pressure 

All the results we describe in what follows were obtained using 
perturbation methods choosing Ti > T\ and using e = {T2 — Ti)/Ti 
as perturbation parameter. We solve the system of equations up to 



e and try to express the nonequilibrium pressure Pyy in terms of the 
density, temperature and the heat flux Qy. Note that here Pyy is con- 
stant. We manage to do this up to third order in e, which corresponds 
to situations not very far from equilibrium. The functional expression 
for the nonequilibrium pressure Pyy turns out to be 



J 



Pyy[n{y),T{y),Qy]=n{y)T{y) 
where Xp is the following constant 



i + 2Xc^n{y) 



n{y)2T{y)^ 



(60) 
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''^"Kn 128 



and where 



A[;i(y),r(y)] = 1 +«(y)(l + 56r(.v));c.— + 



Substituting 5 = in Eq. l l56t yields the functional expression for the 
pressure Pyy, 



Pyy{n,T,Qy) = n{y)T{y) 



1- 



Q 



128 n{y)^T{y)^ 



(63) 



Therefore Eq. i60\ is a natural extension of Eq. J63b . For hard 
spheres the expression which correspond to Eq. l l63b was already 
given in ITlll . For hard disks Eq. l |63t is given for the first time but it 
can be obtained from the appendix in (sj, where a gas of inelastic- 
disks is studied. 



C. Conunent 

The steady state thermodynamics (SST) proposed by Sasa and 
Tasaki t2^ assumes that all nonequilibrium variables and thermo- 
dynamic functions can be expressed in terms of the density, tem- 
perature and the steady state fluxes (heat flux in the case of our ex- 
ample). We found that using Grad's moment method flTl fT^ . this 
condition is satisfied in the simple heat conductive case for the di- 
lute and dense gases which obey Boltzmann's or Enskog's kinetic 
equations in the range of the approximation introduced in Sec.lVlif 
the system is not too far from equilibrium. Unfortunately, because 
kinetic theory defines only the internal energy density, ii{y) 0= T{y), 
and this quantity is interpreted as the temperature, it is impossible, 
within the framework of kinetic theory alone, to construct nonequi- 
librium thermodynamic functions like the Helmholtz energy. Using 



2356 3 3 + 751;rKn2' 



6 

3 + 359;r 



256;^.- 



Kn 



256 



8'- 



32 + 305;r , 3(257;r-3) , 
+784r(v) + '-T{yf 

I 



(61) 



(62) 



Boltzmann's entropy: s = —/In/ {s is the entropy density), one can 
define the entropy and a temperature through (ds/du) = 1/9 where 
9 is the nonequilibrium temperature Bql in what is called extended 
irreversible thermodynamics (FIT) ll3ql . Furthermore one could de- 
fine the nonequilibrium Helmholtz free energy and chemical poten- 
tial, although we have not seen anyone developing such a formalism. 
However within kinetic theory it is not clear which is the difference 
between the temperature and the internal energy density. Even if 
the entropy of the system were given it is not possible to determine 
Helmholtz free energy and therefore it is impossible to determine 
the chemical potential. Using an assumption of SST (2^ that the 
Maxwell relation is satisfied by the nonequilibrium thermodynamic 
variables, then from Eq. l l60b a chemical potential follows without 
having to determine Helmholtz's free energy, although this chem- 
ical potential still has an undetermined part which depends on the 
heat flux and temperature. In a following article IstIi we propose 
a method to get the chemical potential partially using the frame- 
work of SST 12911 . The undetermined part is measured experimen- 
tally. By this method, one can determine the value of the chem- 
ical potential at the interface between equilibrium and nonequilib- 
rium sides. In other words, we will propose an experimental method 
for determining the values of the nonequilibrium variables at the 
equilibrium-nonequilibrium interface following the Sasa and Tasaki 
scheme. Once a chemical potential is given it is possible to com- 
pletely express the condition that the perfect /l wall must have (2^ 
and to discuss about the prediction of (29I1 beyond any doubt. 
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VIII. SUMMARY AND FINAL REMARKS 

In this article we iiave given in detail the hydrodynamic equations 
for a bi-dimensional dense gas of hard disks that follow from En- 
skog's equation using Grad's moment expansion method. Bearing in 
mind that without a strategy to simplify the calculations no manage- 
able hydrodynamics can be obtained we have introduced in Sec. Ivl 
a simplifying scheme that leads to a hydrodynamics which is more 
complete than the one that would follow from a linear approxima- 
tion scheme. As far as we know this is the first time that Grad's 
method used to obtain extended hydrodynmaic equations from En- 
skog's equation has been published beyond the a linear approxima- 
tion in two or three dimensions. 

We have applied this hydrodynamic equations to discuss a simple 
but subtle ID case of heat conduction. The example has the virtue of 
giving the opportunity to discuss the adequacy of the starting point 
of SST 12911 . We have shown that this simple nonquilibrium steady 



state system, not too far from equilibrium, behaves as if nonequi- 
librium thermodynamic variables could be expressed in a functional 
form. Hence, it may be possible to construct a nonequilibrium ther- 
modynamics without any information of the microscopic dynamics 
of the system. First, we have shown this for the dense hard disk 
system. It is well known that for the hydrodynamics obtained from 
Boltzmann's equation using Grad's method this is correct fTTl[T2ll . 
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